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Abstract 

We study the effect of the topology of universe by gauging the 
non-relativistic particle model on the torus and 3-torus, using the 
symplectic formalism of constrained systems and embedding those 
models on extended phase-spaces. Also, we obtain the genera¬ 
tors of the gauge transformations for gauged models. Extracting 
the corresponding Poisson structure of the existed constraints, we 
show the effect of the topology on the canonical structure of the 
phase-spaces of those models and suggest some phenomenology to 
prove the topology of the universe and probable non-commutative 
structure of the space. In addition, we show that the number of 
large extra dimensions in the Phase-spaces of the gauged embeded 
models are exactly two. Moreover, in the classical form, we talk 
over MOND theory in order to study the origin of the terms ap¬ 
peared in the gauged theory, which modify the Newton’s second 
law. 
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1 Introduction 

1.1 Why torus? 

The torus universe model or the doughnut theory of the universe, is the 
model which describes the universe as a doughnut, having surface with 
topology of a two dimensional torus. Historically, the first explanations 
of the shape of the universe were proposed in the mid 60s, after the 
discovery of the CMB by Starobinsky and Zeldovic [T] . 

In experimental point of view, data of cosmos radiation measure¬ 
ments gathered by the satellite CODE, shows small discrepancies in 
temperature fluctuation. This shows that the universe consists of re¬ 
gions of varying densities. Stenemse and Silk proposed that this para¬ 
dox, i.e. the isotropic universe with different regional densities, suggests 
that universe may have a complicated geometric structure [2]. In other 
words, these fluctuations show that multiply connected universes are 
possible, and the simplest and the most probable one is a 3-torus 
Also, simulations of CMB map and the angular power spectrum of the 
temperature fluctuations, considering the torus topology, and compar¬ 
ing them with the observations of the COBE satellite in order to obtain 
the lower limit of universe size, suggest that we live in a small universe 
with the probable topology of the torus laiUElElIH]. 

Another direct observation to detect the topology of the CMB maps 
is an approach which the topology of the universe and the probability 
of being torus has been investigated and called ’’Circles in The Sky” 

piinKiiKiaiis]. 

Moreover, data gathered by WMAP satellite shows more intense 
CMB across one plane of the universe in comparison with others, which 
forms a straight line in the universe. Where radiation surpasses its quota 
for the size of the plane seen, one can say that the universe has over¬ 
flowed in that direction and creates a plane in other directions. Thereby, 
the invisible loop of a torus may have been created perpendicular to the 
direction of the plane. Thus, the analyzed CMB maps from data ob¬ 
tained from WMAP has released some results in favour of a torus form 
of the universe Measurements of WMAP shows that the 

universe is flat with only 0.4% margins of error. On the other hand, flat 
universes with boundaries or edges are not desired mathematically, and 
thus, they are excluded from consideration. Although there are some 
finite compact universe models without boundaries, the torus universe 
is the only one which is both explains a flat and a finite universe |16] . 

Theoretically, string theory and also theories considering extra large 
dimensions suggest that we live in a universe with higher dimensions of 
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space-time and most of the modern cosmological models are founded on 
such assumptions. Moreover, the problems of the standard cosmology 
are avoided, considering higher dimensional spacetime, and also most 
of the predictions of the inflation cosmology are fulfilled via these ap¬ 
proaches [ElIIHlIlg]. 

In order to combine topological theories and extra large dimensions 
universe, it has been shown that cyclic universe models can be acquired 
in a toroidal spacetime which is embedded in a five-dimensional bulk 
with large extra dimensions, and the three dimensional space has been 
shown as a closed ring, moving on the surface of the torus [201121]. 

If we expect that the universe has a topology of a torus, we can con¬ 
struct a gauge theory, using the Lagrangian of a particle on the torus, 
and quantize such a gauge theory, and extract its gauge transformation 
relations. Our goal to study the motion of a non-relativistic particle 
on a torus and gauging that model is to obtain a configuration space 
with extra dimensions. As we know, studying the motion of a free par¬ 
ticle is the most available laboratory, in which we can test whether the 
torus universe has been existed or not. Making such gauge theories and 
studying its Hamiltonian spectrum may help us to understand the real 
topology of the universe. Moreover, with investigating the final obtained 
phase-space, one can check the commutativity and non-commutativity 
of the universe. In addition, we can determine the ratio of two diameters 
of the torus. 

Another point of view, in which we can study the constructed clas¬ 
sical theory on a torus is the MOND theory, that talks over the cor¬ 
rections added to the Newtonian classical mechanics. In the common 
Poisson structure, Hamiltonian equations of motions and Newton laws 
are equivalent |22] . In this article, we construct a classical theory which 
has an unusual Poisson structure due to its constrained structure. This 
Poisson structure adds some additional terms to the Hamiltonian and 
consequently to the equations corresponding Newton’s second law, which 
can be studied via the MOND phenomenological theory. 

Our tool to construct a gauge theory which reduces to a particle on 
the torus after gauge hxing is the symplectic gauge analysis approach, 
will be discussed later. 

1.2 Gauge theories and constraints 

As we know, gauge invariance is one of the most significant and practical 
concepts in theoretical physics. This concept is the cornerstone of the 
standard model of elementary particles. Gauge invariance is due to the 
presence of the important physical variables which are independent of 
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the local reference frames |23j . Whenever a change is applied in an 
arbitrary reference frame, which makes changes in such variables, the 
gauge transformation is occurred. Such physical variables are called 
gauge invariant variables. 

Generally, we deal with gauge invariance, or in the other words, local 
invariance, which produces gauge bosons in fundamental interactions. 
As a physical law, the existence of (local) gauge symmetry in particle 
physics is the sign of the presence of interactions [24j . 

It is very important to know that the quantization of gauge theories 
entails a particular prudence, because of the presence of gauge sym¬ 
metry, which exist some nonphysical degrees of freedom, that must be 
eliminated before and after the quantization is applied [25]. 

On the other hand, in a gauge theory, the equations of motion are 
not able to determine the dynamics of the system thoroughly at every 
moment. Thus, one of the most particular features of a gauge theory is 
the emergence of arbitrary time dependent functions in general solutions 
of the equations of motion. The emergence of such time dependent func¬ 
tions is companied with the relations between phase-space coordinates, 
which are called constraints |26] . 

In order to quantize such systems, the identities between phase-space 
coordinates are classified into two main groups, by Dirac m- The first 
group are identities which are present in the phase-space, similar to a 
coordinate or a momentum variable. These identities, which transform 
the physical system without any changes in the phase-space, are called 
hrst-class constraints, and according to Dirac’s guess are generators of 
the gauge transformations in the phase-space. The second group are 
not related to any degrees of freedom and must be removed. Presence of 
such identities, which are called second-class constraints, indicates the 
absence of the gauge symmetry in the system. Therefore, to gauge a 
system, containing second-class constraints, we must transform them to 
first-class ones, as a first step [28l [29] . 

There are some approaches to perform such a conversion, like BFT 
method [sniiHiiiMiitiiiai], the symplectic formalism [2511M1 IMIEt], 
and the Noether dualization technique [Ml EH SO]. As we mentioned 
before, in order to gauge a system with second-class constraints, we use 
the symplectic approach in order to embed a non-invariant system in an 
extended phase-space nnsasH]. 

1.3 Symplectic Formalism 

Symplectic formalism was introduced by Faddeev and Jackiw [35], to 
avoid consistency problems which spoil the Poisson brackets algebra 
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and consequently fail any quantization techniques in constrained systems 
[441145j . The mathematics of this formalism is based on the symplec- 
tic structure of the phase-space, and therefore, is different from other 
approaches. Also, in the symplectic formalism there is no distinction 
between the first and second-class constraints as in the case of the other 
quantization procedures m- 

The starting point of the symplectic approach is a Lagrangian which 
is first order in the time derivatives. All second order Lagrangian terms 
can be converted to first order ones by enlarging the corresponding con¬ 
figuration space so that it includes the conjugate momentum of the co¬ 
ordinate variables [l6]. Being dependent only on first order Lagrangian 
makes the symplectic approach independent from the classification of 
the constraints into primary, secondary, etc. [39]. In this approach, in¬ 
stead of solving the constraints, one adds their time derivatives to the 
Lagrangian and considers the corresponding Lagrange multipliers as ad¬ 
ditional coordinates m- Also, to convert the nature of second-class 
constraints to first ones, the phase-space would be extended with the 
help of Wess-Zumino variables [50|. After such a conversion, choosing 
conventional zero-modes which are the generators of the gauge trans¬ 
formations and obey particular boundary conditions, one can eliminate 
Wess-Zumino variables, which makes the gauged model equivalent to 
the original system m- 


2 Gauging a non-relativistic particle model on 
the torus 

2.1 Particle on the torus 

In the first part of this article we assume a non-relativistic particle on 
a torus in a three-dimensional configuration space. Considering this 
model, the particle lives on a two-dimensional configuration space, effec¬ 
tively. After our gauging process, we will see that at least one dimension 
is added to the previous configuration space, which makes its space more 
realistic. 

In all sections of this article we get the radiuses of the torus 1 and ?, in 
order to use dimensionless coordinates. Thus, in spherical coordinates, 
the surface of a torus is defined by 

X = (1 ? cos 0) cos (/J 

y = (1? cos 0) sin(/7 

z = 


? sm 


( 1 ) 
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The surface of the torus is described by primary constraint 4>i = 0 
in configuration space for free particle on it. 

(/)i(r, 0) = - 2?cos6'- (1 + ?^). (2) 


In this coordinate canonical Hamiltonian for unit mass is 


He — —{pi + 


pI 


pip 


^ sin^ 9' 


(3) 


In formal constrained analysis we arrive to secondary (finally) constraint 
in Phase-space as 


(j)2{r,6,Pr,P0) = 2 {rpr + 


(,pe sin 9 


(4) 


The set of constraints form a second-class system with non-constant A 
matrix as 


Ai 2 = 4(r^ -F 



(5) 


which makes its embedding by BFT method problematic. This is the 
reason that we use the symplectic approach which is not affected by the 
Poisson structure of second-class constraints. 


2.2 Symplectic analysis of a particle on the torns 

Constructing first-class models from a singular Lagrangian is more straight¬ 
forward in the symplectic formalism than other similar approaches. This 
is done by embedding the primary model in an extended phase-space. 

In this model, the singularity nature of the free particle Lagrangian 
due to its conhguration constraint, (pi{r,9), can be imposed by a new 
dynamical variable (say undetermined Lagrange multiplier) A, in such a 
way that adds the constraints to the free Lagrangian, 

= rpr + 9pe + ifpp, - He - Xi4>i{r, 9). (6) 

Note; In this article the Greek indices, a , f3, a, and (5, are used to 
determine the phase-space variables. 

The symplectic variables and their conjugate momenta as the sym¬ 
plectic one-form can be read off from the model straightforwardly, 

= {r,9,p,pr,pe,Pp,X), 

= iPr,P9,Pp,0,0,0,0). 


( 7 ) 
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The symplectic two-form is defined by 

faf^ = daAf (8) 

Thus, the zeroth-iterated symplectic two-form, using the equation ([SD 
and the symplectic variables and the corresponding conjugate momenta 
(0, will be obtained as follows: 


( Osxs —Isxs OsxlX 

IsxS 03x3 Osxl j . (9) 

0lx3 0ix3 0 / 

This matrix is singular, and so, it has the following null vector, 

nW = (0ix3 0ix3 1). (10) 


Using the zero iterative potential, 

= F, + Ai,/.i, (11) 


the hrst constraint ([2]) will be obtained from following formula. 


4>i=n 


(0) 

a 


dvio) 

Q^{0)a ■ 


( 12 ) 


Substituting the first constraint, obtained from (jl2p into the origi¬ 
nal Lagrangian, we can put the constraint into the kinetic part of the 
Lagrangian. It means that we make the primary constraint 4>i as a mo¬ 
mentum conjugate to the variable Ai. In other words, we convert the 
strongly nonlinear constraint, 4>i, into the momentum (linear constraint) 
of the phase-space. Hence, the first iterative Lagrangian will be obtained 
as 

= rp^ + Ope -h ifp^ - Xi4>i - He- (13) 

We see that the constraint is omitted from the potential. So, for the 
first iterative potential we have. 


= H,. 


(14) 


Now, we read off new symplectic variables and one-form from (I13|) . 

= {r,e,ip,pr,pe,P^,Xi), 

A^a^ = {pr,Pe,P^,0,0,0,4>i). 


( 15 ) 



The corresponding symplectic two-form is constructed as, 


/ 03x3 -l3x3 ufx3\ 

fal3 = I ^3x3 03x3 O 3 XI 1 , (16) 

\ —Uix3 0ix3 1 / 

which; 

Ua = (2r 2?sin0 O) (17) 

The two-form (1161) is a singular one and it has following null vectors, 

= (0ix3 uix3 0 ) , 

4a = (0ix3 0ix3 1) . (18) 

From linear algebra, we know that the linear combination of these null 
vectors is also a null vector. 


na = 4a + /i4a (19) 

Using (jl2p . we obtain the second constraint. 

, , sin 0, , , 

(t)2 = 2{rpr H-^—). (20) 

Now, the second iterative Lagrangian is 

^(2) ^ ^ 21 ) 

and new symplectic variables and one-form are 

= ir,9,ip,pr,pe,P^,Xi,X2), 

^a^ = {Pr,Pe,P.fi,0,0,0,(f)i,(f)2), (22) 


with which we can construct the following symplectic two-form. 


r(2) _ 
J a0 


0 

CO 

X 

CO 

— l3x3 

T 

Ulx3 

vrx3\ 

l3x3 

03x3 

O 3 XI 

wfx3 

-Uix3 

0lx3 

0 

0 

\-Vlx3 

-Wix3 

0 

0 / 


and, V and w are row matrices which are defined as fallow, 


= (2(p, - 2i£^) 22 ^ 0 ), 

Wa = {2r ^ 0). 


(23) 


(24) 
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The corresponding symplectic two-form is non-singular. Thus, it does 
not have any null vector and consequently the iterative process stops 
and no other constraint will be obtained. 

Now, we start the symplectic embedding procedure to convert second- 
class constraints to first ones. The main idea of this procedure is to 
adjoin the Wess-Zumino variables to the original phase-space [50]. In 
order to do that, we expand the original Phase-space by introduction of 
a function G as WZ Lagrangian, depending on the original phase-space 
variables and the WZ variable ct, as the expansion in terms of the WZ 
variables, defined by 


G{r,e,ip,Pr,pg,p^,Xi,a) = (25) 

n=0 

This function is gauging potential and satisfies the following boundary 
condition by vanishing 


G{r,e,ip,pr,pe,p^,Xi,(T = 0) = 0. (26) 

Introducing the new term G into the original symmetrized Lagrangian 
m, we obtain a Lagrangian which depends on both original coordinates 
and WZ variables, 

-I- Lwz, 

= + G{r,e,<f,pr,pg,p^,Xi,a). (27) 

By extending the phase-space, symplectic variables and one-form will be 
extended as 


= ir,0,(p,pr,pe,p^,Xi,a), 

■^5^ = iPr,P9,P<p,0,0,0,4ii,0). 

Galculating corresponding symplectic two-form we have 

7(1) ^ ( fa^ 07xA 

^x7 0 

which has the following zero modes 


II 

(»S 

0 

II 

(41 

») 


(28) 


(29) 


(30) 
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These null vectors are the generators of gauge symmetries, since their 
contraction with the gradient of the potential does not produce any 
constraint [25]. We can use the linear combination of these null vectors, 

na = + hn^2& ( 31 ) 


In order to compute Lwz, we must be assured that no other con¬ 
straint is produced. This mandatory condition generates an iterative 
system of differential equations, defined by the following equation. 


(9^(0)“ da 


(32) 


Substituting (fTTl) into ([32|) . we determine Q^) after an integration pro¬ 
cess as 


e(.) = (2r^ + 


(33) 


Putting into (|27|) . the first iterative Lagrangian will be obtained. 
Hence, for the first-iterated potential we have, 




Using the (j32p for the second time to get we will have. 


(34) 


gi2) 


-2(r2 + 


sm 


-)a2 


(35) 


Substituting into the first iterative Lagrangian, we will obtain the 
second iterative Lagrangian. Consequently, the second-iterated poten¬ 
tial is 

. (36) 

Again, using (1321) to obtain G^^\ we will see that = 0, and so, 

the zero-mode m does not make a new constraint. In conclusion, all 
correction terms G^^\ with n > 3 are vanished. Thus, the gauge invari¬ 
ant canonical Hamiltonian, which had been defined as the symplectic 
potential, is obtained from 

H(c) = +G{r,e,ip,pr,pe,p^,Xi,a), 

= H, + Xi(P,-G^^^ -G^^\ (37) 


and the gauged Lagrangian (1771) . will be 


= L^^'> + G^^'’ + G^^^. 


(38) 
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The generators of infinitesimal gauge transformations can be ob¬ 
tained, using Eicj/, which (/>* are first-class constraints [521 I53j . Also, 
substituting zero-modes (1301) in the following relation. 


= £inf\ 


(39) 


one can obtain the following infinitesimal gauge transformations |44[ I49j . 


6 r = 0, 6 pr = 2rei, 

66 = 0, 6 p 0 = 2?ei sin0, 

dp = 0 , 6p^ = 0 , 

(5A = 62 , 6(7 = Si, 


(40) 


which Ei are infinitesimal time dependent parameters.Obtaining a non¬ 
linear hrst-order Lagrangian, we have constructed a gauge theory with 
the corresponding nonlinear generator functions of gauge transforma¬ 
tions for the model. Thus, the gauge symmetry of the model is deter¬ 
mined via these transformations. In other words, the gained model is 
invariant under these transformations. 

To obtain gauge symmetries of the model, one can use the Poisson 
brackets of the first-class constraints and symplectic variables via the 
following relation [52l [53] , 

= (41) 

Apparently the results obtained from (I41I) is the same as the infinitesimal 
gauge transformations (|40l) . 

Considering constrained analysis of the Lagrangian (|38h and segre¬ 
gating its corresponding constraints in the following section, we study 
the gauge symmetry of the model more easily. 


2.3 Constraint structure of the gauged Lagrangian 


Using the symplectic method, we enhance the gauge symmetry of the 
primary model. In following, we derive constraints and phase-space 
structure of the gauged Lagrangian (l38]l . In this gauged model, new 
dynamical variables A and a appear first-orderly in the Lagrangian. So, 
their momenta are primary constraints in the phase-space. Thus, 


am ^ 

d\W 

am 


0 Pi =px, 

0 /92 = Pa¬ 


so, the total Hamiltonian, corresponding to Lagrangian ([381) . is 


(42) 


ht = He+Pi¬ 


rn 
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In the chain-by-chain method m, the consistency of each individual 
constraint, i.e. pi and p 2 , starts a chain and gives the next element of 
that chain. Also, the consistency of second-class constraints determines 
some of Lagrange multipliers, , while the consistency of Hrst-class ones 
leads to constraints of the next level, 

0 = {pi,HT}, 

0 = {pi,Hc}+ uj^pi,pj}. (44) 

We see that primary constraints are Abelian, i.e. {pi,Pj} = 0. So, 
we arrive to secondary constraints ijji = {/9j,ifc}) where V’l = </*! and 

i>2 = <p2- 

The consistency of second level of constraints gives no new con¬ 
straints, Since, 


= (45) 

The first relation is identically true on the constrained surface, and 
the second one determines a Lagrange multiplier due to the fact that 
{'4’2,P2} / 0. 

All in all, we have the following chain structures. 

Pi ^ V’2 ^ X , 

P2 ^ "02 ^ X . (46) 

Calculating all Poisson brackets, we see that pi is a first-class con¬ 
straint. The Poisson bracket matrix of other four constraints is non¬ 
singular. Non-vanishing elements of that matrix are 


{P2M = ^{^ + Wl), 

(47) 

{V'1,'02} = -^i4: + W2), 

(48) 


in which, tci and W 2 are defined as the components of the row matrix 

(IMl). 

Since the matrix of Poisson brackets is a square matrix with odd 
dimensions, it is a singular matrix. This singularity shows that we have 
more than one hrst-class constraint in our model, other than pi. Re- 
dehning those constraints and requesting hrst-class conditions, we will 
obtain one extra hrst-class constraint, and in conclusion, there will re¬ 
main just two second-class constraints. 
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So, the second first-class constraint is the linear combination of p 2 
and 'i/^i, as <1*3 = /92 + "01 • This first-class constraint strongly commute 
with the two remained second-class constraints, as same as /?i, 

{d'3,V'i} = {d'3,V'2} = 0. (49) 

Therefore, one can rewrite all constraints in the following notation, 

= Pa, 

= V^2, 

^3 = /02 + ^l^^- (50) 

which and $3 are first-class constraints, and and < 1 > 2 ^^ are 
second-class ones. 

Now, we put all second-class constraints into the Hamiltonian to 
calculate the corresponding Dirac brackets. Also, we take first-class 
constraints intact, because they obey the Abelian algebra. 

He = Hc +(51) 

By counting the dimensions of new variables in extended phase- 
space, we find that [Ai] = {Length)~‘^, and [u] = {Length)~‘^. Thus, 
redefining the following variables with length scales, 

Ai = A'“^, (j = (52) 

and replacing them in the Hamiltonian (I5ip . we have, 

ft = Hc+T4.<'l + ^4,l», (53) 

We see that two variables with length dimensions have been added to our 
phase-space. Hence, these length scales extend our configuration space 
from three to five. As a matter of fact, \' and a' can be interpreted as 
large extra dimensions, which are added to spatial part of the phase- 
space, via the potential which carries them in the Hamiltonian. This 
result, i.e. having two large extra dimensions, is in a good accordance 
with [48] . 

2.4 Quantization of the primary model and the gauged 
model 

Taking into the account two primary constraints of the original model, 
(/)! and ( 1 ) 2 , and two second-class constraints, and of the gauged 
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model, and calculating their corresponding Poisson brackets matrix, we 
will have. 




0 

[A + wD 



(54) 


which Wi are the components of the row matrix (|24p . 

In order to determine all Dirac brackets of the original and gauged 
model, we put the inverse of Ajj, in the following formula. 




Non-vanishing Dirac brackets, using components of the row matrices 
(|17p and (I24p . which are common in both primary and gauged models 
are 


Dirac Brackets 

Primary Model 

Gauged Model 

{r,Pr) 

1 

A+w'^ 

1 ^ 

A+w'^ 

{d,Pe) 

^ + ITT? 

1 1 2 ?i )2 


1 

1 

{r,Pe) 

ui(4+wf) 

ui(A+w^) 

{e,pr) 

4:W2 

ui(4+w'-^) 

iiW2 

ui{A+w'4) 

{Pr,Pe) 

4V2 —UIVIW2 

U 2 [Awia{ j—w '2 )—vi+2aV2]—V2 

ui{A+w'j) 

'w 'i , 0 , 

{o', Pa) 

N/A 

1 

(A,Pa) 

N/A 

1 

{Pr,Pa) 

N/A 

2 ui 

{Pd,Pa) 

N/A 

-2 u2 


Table 1: Dirac brackets between extended phase-space variables of the 
torus model 


It is evident that we have obtained some non-commutativity in mo¬ 
mentum parts of the phase-space in quantization process. This non¬ 
commutativity is due to the quantization on the curved space. 

Expanding the Poisson brackets of the gauged model with respect to 
the ratio of radiuses of the torus, and considering the ? —)• 0 , the effect of 
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the topology on the Poisson structure of the phase-space can be studied. 


{r,prY ^ 

{O^PeY ^ 
{r,PeY ' 
{0,PrY '■ 

{Pr,PeY 

{pe,PaY 


1 + 0{Y), 

i + o{Y), 

2? sin 6 


r 

2? sin 


+ o , 

- + 0(e), 


2<j [—pe cos 6 — rpr sin 6 + AYa sin 0) 
—4?sin0 -|- O . 


+ o(.^), 


(56) 


As we see, these Dirac brackets do not have the common canonical struc¬ 
ture. This deformation from canonical structure is due to the change in 
the topology of the configuration space from M” to T^. 

If our test particle can feel such a deviation, we can do some phe¬ 
nomenological research, i.e. we can compare the deviations in the right- 
hand side of the equations (1561) with the discrepancies obtained from 
cosmological data, like the spectrum of energy level of the free particle 
to check whether it is quantized or it is continuous, in order to determine 
the topology of the universe. 

Also, by characterizing hrst-class constraints and Dirac brackets of a 
classical system, its quantized model, say Hilbert space of the quantum 
states, is fully available at tree level, according to Dirac prescription, 

{A, BY ^ ^[A, B], ^Fc\phys>=0, (57) 

which (j)pc is a quantized version of the first-class constraint. Thus, 
due to (j56p in quantized model we derive a non-commutative structure in 
the momentum part of the phase-space. Hence, the quantum mechanics 
obtained from these commutators, is another laboratory to investigate 
the truth of the relations (I56p . and in conclusion, the hypothesis which 
says that we may live on a torus. Also, we can extract some phenomenol¬ 
ogy from these equations to check the non-commutativity structure of 
the space. 


3 Gauging a non-relativistic particle model on 
the 3-torus 

3.1 Constraints of a particle on the 3-torus 

As we have mentioned before, a scenario for the universe as a whole is 
the boundary of a four-dimensional 3-torus. In this section, we consider 
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a test particle on a 3-torus and repeat the previous calculations to derive 
a gauged model for phenomenological purposes. 

The surface of the 3-torus describes the primary constraint = 0 
for a particle living on it. 



(58) 


This 3-torus is described by three dimensionless radiuses as 1, Ci, and 
? 2 - One can reduce this constraint to a polynomial by eliminating the 
radicals (See El). 

The corresponding canonical Hamiltonian for unit mass is 



(59) 


3.2 Symplectic analysis of a particle on the 3-Torus 

Here, like the torus model, we have a polynomial which its consistency 
will determine Lagrange multiplier. Hence, the constrained structure of 
this model is similar to the previous one, and the only difference will be 
occurred in the explicit form of the constraints. Thus, the corresponding 
gauging process for these two models are similar, but because of the 
greater configuration space for 3-torus, we will obtain a new Poisson 
structure. 

The free particle Lagrangian is singular, due to its configuration 
space, affected by a new dynamical variable as an undetermined La¬ 
grange multiplier. A, which adds a constraint to the free Lagrangian. 


l(o) = xpx + ypy + zpz -h sps - He - \i(l)i{x,y,z,s). (60) 


The symplectic variables and symplectic one-form can be read off from 
the Lagrangian, 


= {x,y,Z,S,Pa:,Py,Pz,Ps,>^): 

= {Px,Py,Pz,Ps, 0,0,0,0,0). 


(61) 


Starting the symplectic procedure, the corresponding symplectic two- 
form is obtained as, 



(62) 
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This matrix is singular and so, it has the following null vector, 

n® = (0ix4 0ix4 1) . (63) 

Using the zero iterative potential, 

= 77, + Ai^i, (64) 


the first constraint will be obtained from following formula. 


(t>i = 


ay(o) 

Q^(0)a 


(65) 


which gives the (f58l) . 

In order to remove the constraint from Hamiltonian and add it to 
the kinetic part of the Lagrangian, we substitute (j65p in the Lagrangian 
dsni). As a result, the first iterative Lagrangian will be obtained as, 


= xpx + ypy + zpz + sps - Xi4>i - He, 
and the first iterative potential will be, 

= He. 


( 66 ) 


(67) 


Then, new symplectic variables and one-form are 

= {x,y,Z,S,Pa,,Py,Pz,Ps,X), 

= iPx,Py,Pz,Ps,0,0,0,0,(j)i), ( 68 ) 


which gives the corresponding symplectic two-form. 


and. 


( 04x4 
l4x4 

—77ix4 


— l4x4 
04x4 
0lx4 


04x1 1 ) 

1 / 




d4>i 


(69) 


(70) 


which is the spatial component of the symplectic phase-space, i.e. 
X, y, z, s. Since, the tensor (l6^ is a singular one, it has following null 
vectors. 


hS = (0ix4 Uf, 0) , 
h-2o ~ (0lx4 0ix4 l) ■ 


(71) 
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The linear combination of these null vectors is also a null vector for (j69p . 

fia = hS + hn^2a ■ (72) 

Using (f58]l . we find the second constraint, 

^2 = (73) 

which (9“ is the derivative with respect to the symplectic variables 
Now, the second iterative Lagrangian will be 


L^^'l = xpx + ijPy + zp^ + SPs - Xl4>l - X2(j)2 - He, 

and new symplectic variables and corresponding one-form are 

= (x,y,Z,S,p^,Py,p^,Ps,Al,A2), 

= {px,Py,Pz,Ps,0,0,0,0Al,h)- 

The corresponding symplectic two-form is 


f(2) _ 

1 a/3 


in which, Va and Wa are defined as fallow. 


(74) 


(75) 


04x4 

— l4x4 

^1x4 

Vfx4\ 

l4x4 

X 

O 

04x1 

Wfx4 

—HixA 

0lx4 

0 

0 

—Vlx4 

X 

1 

0 

0 / 


(76) 


V =^. 




(77) 


The two-form (1761) is non-singular. Thus, it does not have any null vector 
and consequently there is no other constraint. 

To start the symplectic embedding process, we expand the original 
Phase-space, using the unknown function depending on the phase-space 
variables and WZ variable, k, which is defined as the following expansion 
with the same boundary condition as 


G{x,y,z,s,px,Py,Pz,Ps,Xi,K) = 


Tl=0 


Introducing the new term G into the Lagrangian 

^( 1 ) ^ ^( 1 ) +G{x,y,z,s,Px,Py,Pz,Ps,Xi,K 


(78) 


(79) 
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extends the symplectic variables as follows, 

= {x,y,Z,S,Pa;,Py,Pz,Ps,Xl,K), 

= iPx,Py,Pz,Ps,0,0,0,0,(})i,0). (80) 


Calculating the corresponding two-form symplectic matrix, we have 


fW ^ ( fad OgxA 

V^1X9 0 ) ’ 


(81) 


which has the following null vectors. 


= 

(»S 

0’ 


4a = 

(42 

o). 

(82) 


These null vectors also generate gauge transformations on the symplectic 
variables ()80p . To continue the procedure, we use the linear combination 
of them, 




= n]z + hhdd 


(83) 


Considering the fact that null vectors (1821) terminates the constraint 
making process, we can use the following differential equation to obtain 


0 


(n) 


Ua 


5y(l) 


(84) 


Substituting (fTTl) into 
variable as, 


9 ^( 0 )“ Ok 

we find as a linear function of WZ 


g(b = n(l,2. (85) 

Putting into (17^ . the potential will be 

= ( 86 ) 

Using the (l8^ for the second time to get g(^\ we will have 

0^'^ =-y Wi>2}. (87) 

Substituting g^^^ into the hrst iterative Lagrangian, we will obtain 

the second iterative Lagrangian with the following potential. 


-g(2). 


(88) 
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Again, using (l8H) to obtain we will see that = 0 and so, 
the zero-mode (|83l) does not make a new constraint. Thus, all correction 
terms with n > 3 are nulls. So, for the canonical Hamiltonian we 
have 

%) = H, + Ai,^i-gW-g(2), (89) 

and for the gauged Lagrangian, 

+g(i)+g(2). (90) 

As we mentioned before, in order to obtain gauge symmetries of the 
model, one can use (HTI) as an option [52l[53], or (l3^ and corresponding 
zero-modes (I82p as another one |44( 149] . which both give the same result. 

5x = 0, 6px = eiUi, 

6y = 0 , 6py = €iU 2 , 

6z = 0, 5pz = eiU 3 , (91) 

5s = 0, 6ps = eiU^, 

5X = €2, 5k = ei- 

These are nontrivial variations that make the system invariant under 
some gauge transformations. So, in the new model, there are some 
generators for gauge transformations which then we seek them. 

3.3 Constraint structure of the gauged Lagrangian of a 
particle on 3-torus 

Now, we can find the constraint structure of the gauged Lagrangian (l90l) . 
using the same method as (I42|l and check consistency conditions. Thus, 
we have 


am 

d\W 

am 

~dkJX) 


0 Pi =PX, 
0 P2 =Pk- 


(92) 


The constraint structure of the 3-torus is similar to the ordinary 
torus. Thus, we have the following chain structures. 


Pl^'tpl^'tp2^X 

P2 ^ i’2 ^ X . 


(93) 


Here, pi is a first-class constraint, while its Poisson bracket and all 
constraints vanish. In order to make another first-order constraint we 
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should redefine them like (j50l) 



^3 = /52 + ^S^\ 


(94) 


which and 63 are first-class, and and $ 2 ^^ are second-class 
constraints. 


Now, we can make the canonical Hamiltonian, 


%) = .H, + 


(95) 


Similar to the torus model, added coordinates to the extended phase- 
space have the following dimensions. 


[A] = {Length) [«^] = {Length) 


Then, we rewrite the gauged canonical Hamiltonian, using the variables 
with length dimension as, 



(96) 


As same as torus model, X' and k' can be interpreted as two large 
extra dimensions which are added to the phase-space. 

We see that the constrained structure for a free particle on the 3-torus 
is similar to the torus. Thus, its Poisson structure and Dirac brackets are 
somehow similar. One can obtain non-vanishing Dirac brackets of this 
model, using (f55]l . Also, we can check the effect of the 3-torus topology 
on Dirac brackets by expanding the Poisson brackets of the gauged model 
with respect to the ratios of radiuses of the torus, considering both 
^ 0 and ?2 0 . 

Like the torus model, we see that these Dirac brackets do not have 
the common canonical structure (See|B]). Thus, we can interpret such 
deformations from canonical structure as the effect of topology on the 
Poisson structure of particle living on a 3-torus. 

3.4 How to extract phenomenology related to MOND 


theory 


As we know, to hnd Poisson structure between two functions A and B 
in a phase-space which owns itself a Poisson structure as )?/?}*) one 
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can use the following relation 




dA dB, 
d(a d^0 ' 


(97) 


which {^ 0 ,^/ 3 }* is the Dirac bracket between phase-space variables iBl 
The relation (f^lTD is used to obtain the time evolution of phase-space 
variables and consequently the Hamilton equations of motion. 

In order to relate the obtained gauge theory on the 3-torus and 
MOND theory [56l [571 ESI EH], we rewrite the Hamilton equations of 
motion for dynamical variables of our gauged model and combine them 
to obtain the Newton’s second law, and more interestingly its corre¬ 
sponding corrections. 

To start with, we consider the general radial potential V{r), added 
to the Hamiltonian which does not change the derived gauged theory. 


jj Pi'P , .rf S 
H = — \-V{r), 


(98) 


which r is the radial vector, i.e. = x‘^ + y'^ + z'^ + , in the Cartesian 

coordinates. 

Hence, due to the nontrivial Poisson structure, the Newton’s second 
law in the direction of xv will be 


Xi' = - 


dV 

dxji 




(99) 


where, can be interpreted as the modification term of the 

Newton’s second law. 

Let‘s start studying this modification in details. As we mentioned 
before, due to the amendments which are imposed by Dirac brackets 
ascribed to second-class constraints, the equations of motions shall surely 
include some corrections. Thus, for the equations of motions we have. 


Xj/ {Xi^,Hc} = {Xi>,^a} 


dH 

Wa 


Xi' ~ {{xii,Hc},Hc}. ( 100 ) 

Calculating the explicit relation of acceleration, i.e. x*/, we arrive to 

977c. rr i977c c . cr. rr. 

Xi' « {Xii,Pi\{^^,Hc\ -7 -^\.{{Xi' ,Pj}-^,^a\ - {{U,Hc\,Xi'\). 

( 101 ) 

It is easy to investigate that the equations of motion in the con¬ 
structed theory give the following correction for the acceleration of the 
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free particle which does not feel any potential. As we know, these cor¬ 
rections are equal to zero for common situations. 



It is evident that these accelerations depend on the position, momentum 
and the kinetic energy of the particle and are of the order of {Length)~^. 

For a particle with unit mass which is affected by the potential in a 
theory with flat space and the topology all the terms in (1991) vanish 
but the first term — So, we obtain the Newton’s second law as 
J" = Xj/. 

On the other hand, for a gauged theory, the relation (jin2jl contains 
extra terms which have the capability to be explained as the correction 
of the Newton’s second law. These correction terms can be used for 
phenomenology of the MOND theory, and somehow, be the candidate 
to explain dark matter. 

As a matter of fact, these corrections are classified into two cate¬ 
gories. First, terms which are added to (|100l) in a proper gauge, to 
convert weak equalities to strong ones, or in other words, changing He 
to Ht, in order to gain full dynamics of the particle. Moreover, for the 
case of having the common Poisson structure, the only survived factor 
which is derived from the first term of (I99p . is —^ 7 - But, in the model 
of the particle, living on the torus, as we saw in the last section, we 
encountered some deviations in the common Poisson structure of the 
phase-space. Thus, aside from the survived first term which adds some 
corrections itself because of the deformed Poisson structure of the phase- 
space, there are also correction terms, arised from other terms. All these 
modification factors can be used to study MOND. 

In addition, if we consider the gravitational potential of the mass M 
as 1 / = then we can find the correction which is imposed on the 
universal law of gravity via the attendance of the particle on the torus 
in tree-dimensional space, and in the presence of the extra gauged coor¬ 
dinates [601 EH [62]. Thus, one can conclude that for the particle which 
is affected by the gravitational potential, despite the added corrections 
for accelerations, other corrections are emerged for the gravitational po¬ 
tential itself. 

Moreover, these modification factors can be interpreted as the func¬ 
tions of the new coordinates X' and H. Prom this point of view, we can 
elucidate these terms as large extra dimensions in Randull-Sundrum sce¬ 


nario 
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Conclusion 

Based on the observations and the theories in cosmology, we assume that 
our universe has a topology of a torus. We consider a non-relativistic 
particle as a test object in this background and find the corresponding 
gauged phase-space in classical mechanics framework to obtain classical 
equations of the particle. Also, using Dirac’s approach, we gain the par¬ 
ticle’s quantum mechanics to the first order of fi and associated Hilbert 
space. 

First, we construct a gauge theory on a torus in three-dimensional 
space as a toy model. We try to concede two degrees of freedom to 
this particle which lives on a two-dimensional world, using symplectic 
embedding approach. We show that the interaction which is added 
to the Hamiltonian of the particle via this process is as the inverse of 
the length to the power of four on the axis of the gauged degrees of 
freedom. We also show that these two added degrees of freedom is due 
to the imposing the constraint of being on the torus on the Hamiltonian 
with the help of Lagrange multipliers, and the gauging WZ variables. 
Extracting its Poisson structure, we can quantize the model, by replacing 
Dirac brackets with quantum commutators. Afterwards, we study the 
effect of the topology on Dirac brackets by enlarging one radius of the 
torus in comparison with another. We witness that the model of a free 
particle in the torus converts to a particle on an infinite ring with the 
topology of 5^, which is the very effect of the topology on the Poisson 
structure of the phase-space. 

With the help of this toy model, in the main part of this article, 
we constrain the particle on the hyper-surface of a three-torus in a 
four dimensional configuration space. Because of the second-class con¬ 
straints which are imposed on such a particle, the particle will live on 
a six-dimensional phase-space or three-dimensional configuration space. 
Then, using the symplectic formalism and embedding the phase-space 
to an extended one, to the particle’s motion degrees of freedom we add 
two more ones. These degrees of freedom are equivalent to two second- 
class constraints, gauged in the hnal model. By constructing the main 
brackets of the phase-space, we obtain classic and quantum mechanics 
of the particle, from which we can interpret deviations occurred from 
classic and quantum mechanics in space due to the topology of the 
torus. Hence, we can say that the effect of the topology 5^ x x 
on Dirac brackets, leads to the presence of the particle on a three di¬ 
mensional manifold, i.e. the particle travels on the surface of the 3-torus 
and is not floated in its bulk. The obtained quantum mechanics can be 
investigated by studying the energy spectrum of the free particle in the 
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model. In the gauged particle’s quantum mechanics which is affected by 
the space topology, main commutation relations change basically. So, 
even in the free particle’s Schrodinger equation of motion, there will be 
an added potential which makes the corresponding energy spectrum dis¬ 
crete. Hence, by investigating the free particle’s energy spectrum and 
comparing with the very spectrum which is obtained via quantum me¬ 
chanics of the model, we can determine the corresponding parameters 
accurately. 

Another research work in cosmology can be done with the help of 
above model to obtain the proposed corrections from MOND theory. 
Using the equations of motion obtained from the gauged Hamiltonian, 
we gain the modified form of the Newton’s second law. The obtained 
correction terms in our model can be compared with the cosmological 
observations for investigating MOND, and at the same time, it verifies 
the MOND theoretically. 

Moreover, the obtained model has two extra gauged degrees of free¬ 
dom comparing to the common three-dimensional space’s degrees of free¬ 
dom, which add interactions proportion to inverse length square and in¬ 
verse of the length to the power of four to Hamiltonian. Particularly, if 
we consider that the test particle is affected by the gravitational poten¬ 
tial of another one in the origin of the space, then in the gauged model, 
the particle faces the corrections due to the two extra dimensions, which 
can be related to the brane world cosmology. 
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A 


As we mentioned before, we reduce the primary constraint of a particle 
model on a 3-torus to a polynomial by eliminating the radicals as follow, 


(^1 = + + Qx^s‘^ 

+ + 62:^5^ + — 12x^s^ + 12x^y^s'^ — 12y^s^ + 12x^2:^s'^ 

+ 12y2^2g4 _ ^222^4 _ 43.2^-2^4 _ 12y2^2g4 _ 12^2<^2^4 ^ ^-2^4 _ 43.2-2^4 

- 4y^?|s^ - 12z^?|s'‘ - + 4x®s^ + Ay^s"^ + 4z®s^ 

- 4?fs2 - 4 ?|s 2 - 12x^5^ -p 12x2y4s2 - Uy'^s'^ + Ux^z'^s'^ + 12y‘^z'^s^ - 12z‘^s^ 

- 4x^?f + 12y‘^qfs‘^ + 12z‘^qfs‘^ + Aqfs"^ - 4x^?|s^ - 4y^?|s^ + 122;^ 

- 60?f ?|s2 + 4?|s^ + 12x^5^ + Ux'^y'^s'^ - 2Ax‘^y‘^s^ + I2y^s^ + I2x^z‘^s‘^ 

+ I2y z s — 24x z s + 24x y z s — 2Ay z s -h I2z s + 4x qs — 12y 

io4-'2 2 o2-2 2 o 22-2 2 |o 2-2 2 © 22-^2 2 o/i 22-2 2 

— I2z c^iS —8x<^iS —8xy<^iS +8y qs — 8x z qs — 24y 2 
+ 8z q5 +4^5 +4a: ^2-^ +^y ^ 2 ^ — I2z <^ 2 ^ —oO<;i<; 2 S — 8x ^ 2 ^ 

+ 8x^y^?|s^ — 8y^<f|s^ — 8x^z^<f|s^ — 8y^z^^|s^ + 8z^?|s^ — 24x^?i?|s^ 

- - 120^;^?i^|s^ + 24?^ + 4<f|s^ - 4s^ - 4x® - Ay^Az^ + 4x^y® 

+ 4x^^;® + Ay‘^z^ - +Ax^ql - 4y^?® - Az^qf - 4?® + 4x^?|’ + 4y^<f| - Az^q^ 

- 4<f<| — 4<| + 6x^ + 6x^y^ — 12x^y^ + 6y^ + 6x^2:'^ + — 12x^2^ 

+ 12x^y^2'^ — 12y^z‘^ + Qz'^ + 6x^qf + Gy'^'ff + Gz'^^f — 12x'^^f — 4x^y^?f 
+ Ay'^qf - 4x^2^+ 12y‘^z‘^^f + Az‘^^f + + 6x^?| + 6y^?| + Gz"^?! + 

- 12x2^1 + I2x‘^y^4 - 1V4 - ^x‘^z‘^4 - 

- 12y\f + 4?i ?| + 6?| - 4x2 + 4x®y2 _ I2x^y2 ^ ^2x2y2 _ 4^2 

+ 4x^2^ + 4^622 _ 12x^2 -p 12x2/22 _ 12// + 12x2^2 12x^y2^2 

- 24x2y2/ _g 12// _ 4/ + 4x®/ _ 4^6-2 _ 4^6-2 _ i2x^/ - 4x2// 

+ 4y\f _ 4x2^4/ _ 12y2^4-2 ^ 4^4-2 ^ ^2x2/ + 4x^y2^-2 _ g^2y2-2 

+ 4y2 / + 4x"^z2/ - 12y"‘2;2^2 _ g^' 2 ^ 2~2 _ g^2y2^2-2 8y2^2-2 4^2-2 _ 4^2 

+ 4x®/ + 4// - 4z®/ - 4// - 12x^/ + 12x2// - 12// - 4x2^;^/ 

- 4y2// + 4// - 4x2// + 12/// - 60/// + 4// + 12x2/ 

+ 12x^y2/ _ 24x2y2/ + 12y2/ + 4x'^// + 4/^2/ - 8x2^2/ 8x2y2// 

- 8/^2/ + 4^2/ + 4x\f / - 12/// - 60/// - 8x2// _ 8x2/// 

+ 8y2// - 24x2^2// - 72y^z‘^^f4 + 24^2// + 4// - 4/ + / + / + 1 
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B 

Dirac brackets of the phase-space variables are shown in the following 
table. As it is evident, these brackets do not include any correction to 
the first order of q, i.e. — Jap + O +0 (^l) +0 (?i? 2 )- Here, 

we dehne + s^. 


{,r 

Qj 

A 

K 

Pj 

PX 

Pk 


0 

0 

0 


0 

0 

A 


0 

0 

0 

1 

0 

K 



0 

0 

0 

1 

Pi 




r'^ 

0 

8gj(r2 - 1)3 

Px 





0 

0 

Pk 






0 


Table 2: Dirac brackets between extended phase-space variables of the 
3-torus model 
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